Supersolid behavior of nonlinear light 
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We present a formal demonstration that light can simultaneously exhibit a superfluid behavior 
and spatial long-range order when propagating in a photonic crystal with self-focussing nonlinearity. 
In this way, light presents the distinguishing features of matter in a "supersolid" phase. We show 
that this supersolid phase provides the stability conditions for nonlinear Bloch waves and, at the 
same time, permits the existence of topological solitons or defects for the envelope of these waves. 
We use a condensed matter analysis instead of a standard nonlinear optics approach and provide 
numerical evidence of these theoretical findings. 



Matter can undergo a phase transition at ultracold 
temperatures in which exhibits a superfluid behavior and, 
at the same time, present all the characteristics of a crys- 
talline solid. This new phase of matter is known as 
"supersolid" and it has been experimentally proven in 
helium-4 in recent years [1 j. Apparently, light is unre- 
lated to these phases characteristic of condensed matter. 
However, analogies between condensed matter and opti- 
cal systems are increasingly appearing in the literature 
[2]. Even the concept of liquid phase of light has been 
already suggested [3 . In this letter, we will take this 
analogy a step further and we will demonstrate a formal 
equivalence between regular light structures propagating 
in a nonlinear photonic crystal and ultracold matter in a 
"supersolid" phase. For the purpose of establishing this 
equivalence, we will use the formalism of condensed mat- 
ter and particle physics [4] instead of using a standard 
nonlinear optics approach. 

In Hamiltonian formalism the nonlinear Schrodinger 
equation for a periodic medium with Kerr nonlinear- 
ity can be obtained from the Hamiltonian density as 
idab/dz = dH/d^ = (-V, 2 + V(x)-#|(/>| 2 )^. In our 
case Vt is the 2D transverse gradient operator and the 
potential V(x) = — (n 2 (x) — n 2 ,) represents the periodic 
modulation of the square of the refractive index with re- 
spect to the reference value n 2 ,. Transverse and axial 
coordinates are normalized. The energy of a propagat- 
ing solution is given by E = j d 2 xH((j)). Now, instead of 
using E we introduce the optical equivalent of the free 
energy in statistical physics F = E — \iP = J d 2 xJ r {qb) 
where \i is the propagation constant of a stationary so- 
lution and P is the system power (whose, in absence 
of losses, is a constant in the propagation). Thus P 
plays the role of N, the particle number in statistical 
physics, and \i plays the role of the chemical potential. 
In terms of the free energy density the equation of motion 
is id$/dz = dT/d^f = (-V? + V(x) - /i - g\4>\ 2 ) 0. It 
is easy to check that the relation between the solutions of 
both equations is simply <fi = (j)e~ lfiz . Thus an stationary 
solution (j) with propagation constant \i is equivalent to a 
^-independent solution that is an extremum of T veri- 



fying dF/dcjf = 0. For the case under consideration, the 
optical free energy would be given by (we use <p instead 
of $): 



F = d 2 x 



Vt<£*V t <£ + y(x) 



(i) 

In this context the typical nonlinear optics P(fi) curve for 
a stationay solution appears as the conventional equation 
of state of statistical physics P = —dF/dfi. 

On the other hand, we are interested in the analysis 
of stationary solutions whose amplitude is invariant un- 
der finite translations of value a, where a is the period 
of the potential. If <p so \ is a solution that satisfies the 
equation for stationary states and that simultaneously 
verifies the condition | (j) so \ (x + a) | = |0 so i(x)| then the 
total potential V so \(x) = ^( x ) — # I ^soi ( x ) 1 2 occurring in 
such equation will be periodic with period given by a, 
V so \(-k + a) = 14oi( x )- Self-consistency implies that so i 
must be equal to one of the Bloch functions so i(x) = 
/q ; /3 ( x ) = x ^Q;/3o ( x ) forming the spectrum of the 
nonlinear operator generated by itself and it will be char- 
acterized by a propagation constant /i = /iQ ; /3 . Using the 
nonlinear operator generated by <p so \ (including the total 
potential V^ Q i( x ) = ^( x ) — ^|0soi( x )| 2 ) we construct a 
nonlinear Wannier function basis by means of standard 
techniques. We assume that all functions are defined 
within a periodic region called the basic domain, of 
spatial dimensions Na x Na so that we represent any 
arbitrary field amplitude at a given axial position z as 
0(x, z) = J2i, a c i,a( z )W£ ( x - x^).The nonlinear Bloch 
solution so i(x) is represented by the site-independent 
and ^-independent coefficients d° l = ^ M e* Q x * where 

77 M = y/P/N. Note that Q is discretized because <fi is de- 
fined in the periodic domain Q so that Q m = 27rm/(7Va) 
where m G Z 2 and fulfills -N/2 < m Xl m y < N/2. We 
introduce now the Wannier expansion into the expression 
for the optical free energy and integrate out the tranverse 
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coordinates x to obtain 
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f = -%)EE4( c ^ 



riv 1 i—rij. 



ql 2 -^£M 4 + (h.o.t.), (2) 



where ni and n2 are the "horizontal" and "vertical" lat- 
tice vectors, t(fi) = -Z^.o+ni = _L o,6+n 2 is the ef- 
fective nearest neighbors coupling and l(fi) = Lqq and 
U(fi) = 2Tqqqq are the two on-site couplings. The second 
and fourth order couplings L^(/i) and XU^(/i) are ob- 
tained as overlapping integrals of the nonlinear Wannier 
basis associated to the stationary solution and, for this 
reason, they depend on (/i, Q,/?o) (f° r simplicity, we dis- 
regard interband interaction terms and elliminate band 
indices): L^) = J n d 2 xWf (-V? + V(x) - M ) 

and T^-(/z) ee ±f n d 2 xgWrWfW[W?. The expres- 
sion (h.o.t.) in ([2| stands for higher-order terms involving 
interactions at longer distances. 

The main difference between ^ and previous ap- 
proaches is that the nonlinear Wannier basis provide dif- 
ferent coefficients than those obtained using linear Wan- 
nier functions or localized single potential solutions as in 
the tight binding aproximation. Since our aim is to anal- 
ize the stability behavior of the nonlinear Bloch solution 
characterized by \i the election of the nonlinear Wannier 
basis associated to it is a natural choice. With this pur- 
pose in mind let us write the z-dependent Wannier co- 
efficients as a(z) = e iQ - x ^i(z) = e iQ,x ? (r^ + A$i(z)) 
to formalize the fact that we want to analyze the dy- 
namics associated to the stationary solution described 



by c\ 



,sol 



Q x ;. In this way the dynamic information 



is encoded in the envelope coefficients &i(z). 

In this Letter we are interested in the regime where 
perturbations have a correlation length larger than the 
potential period — £ ^> a. Physical phenomena in this 
regime present a collective spatial character that will re- 
flected in the properties of the discrete envelope function 
&i(z). Spatial collective effects will be represented by 
values of <^ that will fluctuate smoothly in space. It 
is then natural to take the continuous limit of the dis- 
crete optical free energy by considering the limit a/( — > 
and by introducing the continuous spatial envelope func- 
tion $(xt,z) defined as Q>(pc^z) = <&i(z). We substitute 
Ci = e^'^&i into ([2| and then take the continuous limit 
using standard techniques. In order to simplify the re- 
sult, we consider only "diagonal" nonlinear Bloch waves 
with Q\ = Q<2 = 27rm/(7Va), m G Z. The phase dif- 
ference between two neighboring sites of the nonlinear 
Bloch wave is then A(j) m = 2i\mjN . In order to relate 
the optical effects driven by F to condensed matter and 
particle physics phenomena we need to re-define the op- 
tical free energy in some cases. For this reason, we in- 
troduce the sign-changed optical free energy defined as 




<b) (c) 
Figure 1: Effective potential for different configurations 

F = sign(6 m )F 

d 2 x [|MaO|V$*V$ + zsign(6 m )(v m (/i) • V$)$* 



+^(|<D| 2 ) + (h.o.t.)], 



(3) 



where 6 m (/i) cos(A0 m )t(/i)a 2 , v m (/i) v m (l, 1) 
with v rn (p) = — 2 sin(A0 m )t(/i)a and M m (/i) = 
Z(/i) — 4t(/i) cos(A0 m ). In this case (h.o.t.) in- 
cludes terms with higher-order derivatives and hi = 
sign(6 m ) (M m (fj) |$| 2 - ^|$| 4 ). In principle, the use 

of F or F is irrelevant as far as dynamics is concerned 
since dynamics remains the same under the change F — > 
—F. The advantage in the use of F is that it contains 
in all cases a positive definite kinetic energy term, as 
required in quantum field theory to define a proper vac- 
uum state [5]. By doing this, equivalences are straight- 
forward to achieve even with negative values of b m . By 
means of F we can proceed to determine the nature of 
the ground state of the system by analyzing the behav- 
ior of the so-called effective potential U(<&) "a la Lan- 
dau". Landau theory is a mean-field approach used to 
characterize phase transitions in condensed matter and 
particle physics [4] . The qualitative form of the effective 
potential is represented in Fig. ([I]) for different signs of 
coefficients. One can easily recognize in this figure that 
there are only two configurations for which a spatially 
uniform envelope solution 3>(x) = $o is allowed (cases 
(b) and (d)). They correspond to extrema of the opti- 
cal free energy density dT /d<&* = dU/d<&* = given by 
the condition |3>o| = M m (/i) /U{n). The ground state of 
the system is the state that minimizes the free energy, 
thus, only the two cases in (d) can provide an spatially 
homogeneous ground state. The former analysis is iden- 
tical to that performed in condensed matter physics to 
establish the nature of phase transitions in the mean-field 
regime, here \i playing the role of temperature. In this 
way, the signs of the 6 m (/i), M m (/i) and U(ja) coefficients 
determine the nature of the ground state and, therefore, 
in which "phase" light is. It is important to strees that 
the ii dependence in these "Landau coefficients" is the re- 
sult of using nonlinear Wannier functions. The ground 
state in Fig[TJd) is degenerate since all solutions of the 
type = |3>o| e are minima of the effective poten- 
tial and thus they all have the same free energy. The 
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optical free energy is invariant under a U(l) phase trans- 
formation — > e m <l>. The ground state $0, however, is 
not since this phase transformation maps it into a differ- 
ent degenerate solution with the same free energy. 
This mechanism is well-known in condensed matter and 
particle physics and it is known as spontaneous symme- 
try breaking (SSB). Superfluidity, superconductivity, the 
Higgs mechanism or the chiral phase transition in quan- 
tum chromodynamics are physical phenomena related to 
the same mechanism. The SSB mechanism have distinc- 
tive properties: (i) appearence of a non-zero order pa- 
rameter in the broken phase, (ii) existence of topological 
solitons or defects, (iii) presence of masless or long-range 
excitations (Goldstone bosons). In this language, our op- 
tical system in the Fig[TJd) configuration is in a broken 
phase, a phase in which U(l) symmetry has been spon- 
taneously broken. We expect then to find the optical 
counterparts of the aforementioned properties. In this 
Letter we will pay attention to properties (i) and (ii), 
leaving the analysis of (iii) for a further publication. 

We can now establish a link between the stability of 
nonlinear Bloch waves and the SSB mechanism. A non- 
linear Bloch wave is characterized by an homogeneous 
envelope. According to our previous analysis, only in 
the configurations in Figjljd) this envelope function cor- 
responds to the ground state of the optical free energy. 
On the other hand, in this configuration the system is 
in the U(l) broken phase. Since the ground state is a 
stable solution, i.e., fluctuations around it cannot trans- 
form this solution into a different one by any dynamical 
mechanism, we infer that the stability condition of a non- 
linear Bloch wave is achieved in the broken phase, i.e., 
when (b m > 0, M m < 0, U < 0) or (b m < 0, M m > 0, 
U > 0). The properties of a stable nonlinear Bloch wave 
are then identical to that of a superfluid as far as its enve- 
lope is concerned. However, the solution simultaneously 
presents spatial long-range order. Using quantum field 
theory notation, if |0) stands for the ground state given 
by the nonlinear Bloch wave in the mean-field regime 
at a given fi: (i) T a |0) = e iQ a |0) where T a is the lat- 
tice translation operator and (ii) (0| 4> |0) = |$o|e^ ^ 0. 
Physically speaking, the "soliton lattice" filling the pe- 
riodic medium completely and described by the nonlin- 
ear Bloch wave presents perfect spatial long-range order 
as in a cristallyne solid and, at the same time, it pos- 
sesses a non- vanishing order parameter indicating it is in 
a superfluid phase. The dynamics of low energy fluctua- 
tions around such a solution is determined by Eq. ^ and 
has to show identical features than a superfluid. In this 
way, we demonstrate that, under the specified conditions, 
light fulfills the definition of a supersolid, that is, it is a 
spatially ordered system (like in a solid or crystal) with 
superfluid properties. 

We have performed a number of numerical experiments 
to check both qualitative and quantitative the validity of 
our condensed matter approach. We have modelled a 
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Figure 2: Landau coefficients (6 m , M m , U) in terms of the 
propagation constant \i for unstaggered (dashed line) and 
staggered (solid line) nonlinear Bloch solutions. 
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Figure 3: St ability /unst ability patterns of phase under prop- 
agation: (a) for a perturbed unstaggered solution phase disor- 
ders quickly, (b) for a staggered solution spatial order remains. 



nonlinear photonic crystal with C4 symmetry formed by 
a nonlinear material with refractive index n\ with em- 
bedded circular inclusions of a linear material with index 
n2 < ni, the index contrast being given by the potential 
Vo = — (??| — Tii) > 0. This type of nonlinear photonic 
crystal can be achieved in standard or in chalcogenide 
photonic crystal fibers [6 or in laser-writing waveguides 
[7] . The nonlinearity is self- focussing and of the Kerr type 
Vnl = — 2 (9 > 0). First of all, we have evaluated the 
dependence of the "Landau coefficients" (6 m , M m , U) on 
\i . We have proceeded as follows: (i) we numerically 
evaluate a nonlinear Bloch solution (j) so \ at a given /i, (ii) 
we determine the nonlinear operator associated to (f) so \ 
(i.e., including the nonlinear potential Fnl = — #|0soi| 2 ) 
which numerically will be a matrix, (iii) we find the spec- 
trum of this matrix formed by Bloch eigemodes since the 



(a) (b) (c) 




Figure 4: Supersolid light topological defect for fi — —0.1, 
Vo = 2, and charge +1: (a) simultaneous representation of 
the amplitudes of the envelope <I> and full functions for this 
solution, (b) representation of its phase and (c) P{fj) diagram 
calculated using the envelope equation (solid line) and the full 
equation (triangles). 
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total potential, including the linear and nonlinear part, is 
periodic, (iv) we construct the nonlinear Wannier basis 
out of these Bloch modes using the standard technique 
[8] and (v) we evaluate the "Landau coefficients" using 
their definitions as overlapping integrals of Wannier func- 
tions. A given nonlinear Bloch wave characterized by \i 
then univocally provides a specific value for (b m , M m , 
U). In Fig. [2] we represent the functional form of these 
coefficients in terms of \i for two solutions with differ- 
ent pseudo-momentum Q: a solution with Q = with 
a phase difference between neighboring sites A</> m = 
(unstaggered) and a solution with Q = (ir/a, ir/a) with 
A0 m = 7r (staggered). By analyzing the signs of their 
Landau coefficients we immediately recognize these two 
configurations correspond to the cases (b) (unstaggered) 
and (d) (staggered) in Figfl] The main difference be- 
tween these two configurations arise from the different 
sign of b m , which has its origin in their different phase 
structure since b m ~ t(fi) cos(A0 m ). In the case ana- 
lyzed here t(ja) > for all \i in both configurations, so 
that the different form of the effective potential in both 
cases is a pure effect of the underlying phase of the non- 
linear Bloch wave. According to our condensed matter 
analysis, the unstaggered solution can exist as an homo- 
geneous solution but it corresponds to a metastable state 
that eventually will decay into a different state. It cannot 
be stable. On the contrary, the staggered solution is the 
ground state of a light supersolid and, consequently, it is 
necessary stable. We have performed a numerical stabil- 
ity analysis of these two configurations that confirm these 
predictions. A small initial perturbation of the unstag- 
gered solution gives rise, after a short propagation, to the 
breaking of spatial long-range order crearly reflected in 
the progressive disordering of the phase of the solution — 
see Figj3ja). The staggered solution is, however, immune 
to perturbations and preserve spatial long-range order in 
phase and amplitude — see Figj3jb). Similar examples 
can be found in the literature that can be explained by 
the same mechanism, e.g., stability of large truncated 
nonlinear Bloch waves [9]. As mentioned before, one of 
the most paradigmatic properties of a superfluid is the ex- 
istence of topological solitons or defects. They are holes 
in the superfluid around which the superfluid flows with a 
quantized circulation. Therefore, it is expected that the 
envelope function <E> supports the existence of the optical 
counterparts of these objects in the broken phase. We 
have proven indeed the existence of these type of solu- 
tions by numerically solving the same nonlinear photonic 
crystal structure using: (i) the effective equation for 
associated to F — Eq.([3| — using the values of (6 m , M m , 
U) given in Fig.Q and (ii) the full equation for the orig- 
inal field <j) including the periodic potential V(x). It is 
remarkable that both approaches are in excellent agree- 
ment with each other both qualitatively and quantita- 
tively. As correctly predicted by theory, the topological 
vortex only exists on top of the staggered solution — see 



FigQb) — which corresponds to the ground state in the 
broken superfluid phase |$q| 7^ 0- From an optical point 
of view, one could consider this object as a type of dark 
soliton, however, the striking property of this dark soli- 
ton is that exists in a self-focusing medium when they 
are usually associated to defocusing media. In Figj4ja) 
we can appreciate the excellent fit of the envelope func- 
tion to the solution of the full equation. This excellent 
quantitative agreement is even more clearly appreciated 
in the calculation of the P(fi) curve for a family of topo- 
logical vortices using both approaches shown in FigQc). 
In summary, FigsQa) and (b) show the simultaneous 
presence of superfluid behavior (nontrivial amplitude and 
phase of the topological defect) and spatial long-range or- 
der (perfect staggered order of the background in all the 
domain). Besides, numerical stability analysis indicates 
that these solutions are stable under propagation during 
long distances. These features represent a clear numerical 
evidence of the supersolid behavior of light in nonlinear 
photonic crystals predicted by theory. In this sense, it 
has been shown that light can be treated analogously as 
matter not only in their qualitative aspects but by using 
a condensed matter formalism. We call this approach to 
nonlinear optics photonic condensed matter. 
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